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molten polymer flow is limited by solid surfaces (e.g. in extruders 
or in compression molding). The rheological behavior of a poly¬ 
mer melt may be well known, but unexpected phenomena such as 
sharkskin instability and melt fracture, which are intimately related 
to polymer-substrate interactions, occur at the interface of the melt 
and the mold or extruder.®’^ Polymer chains can adsorb on sur¬ 
faces and modify their effective interaction as well as resistance 
to their relative motion, a phenomenon with consequences in areas 
which range from tribology to biology. It has been theoretically 
proposed that, under certain conditions, a highly entangled poly¬ 
meric fluid may not maintain the usual stick boundary condition 
at a solid wall.^ A significant and unambiguous violation of the 
law of the stick hydrodynamic condition has been experimentally 
observed under stress.®’*® Adhesion or friction occurs and the con¬ 
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A method is formulated, based on combining self-consistent field 
theory with dynamically corrected transition state theory, for estimat¬ 
ing the rates of adsorption and desorption of end-constrained chains 
(e.g. by crosslinks or entanglements) from a polymer melt onto a 
solid substrate. This approach is tested on a polyethyiene/graphite 
system, where the whole methodology is parametrized by atomisti- 
cally detailed molecular simulations. For short-chain melts, which 
can still be addressed by molecular dynamics simulations with rea¬ 
sonable computational resources, the seif-consistent field approach 
gives predictions of the adsorption and desorption rate constants 
which are gratifyingly close to molecular dynamics estimates. 


Introduction 

Understanding the fascinating and complex dynamics of melts of 
large flexible polymer coils close to solid substrates has been an 
ongoing challenge for many decades. From the point of view of 
molecular simulations, the spectrum of length and time scales as¬ 
sociated with polymer melts of long chains poses a formidable 
challenge to studying the long-time dynamics. *’^ The topologi¬ 
cal constraints arising from chain connectivity and uncrossability 
(entanglements) dominate intermediate and long-time relaxation^ 
and transport phenomena when polymers become sufficiently long. 
Atomistic molecular simulations of dense phases of soft matter 
prove to be difficult for many systems across length and time scales 
of practical interest. Even coarse-grained particle-based simulation 
methods may not be applicable due to the lack of faithful descrip¬ 
tions of polymer-polymer and polymer-surface interactions. While 
the study of polymer adsorption from solution onto a solid surface 
has a long history, the adsorption from a melt has been less stud¬ 
ied. 

The thermophysical and rheological properties of polymer melts 
are of paramount importance in plastics processing. Many ther¬ 
moplastic products are formed by means of processes in which the 


formations of macromolecules near the interface are often different 
from the bulk. After cooling and solidifying this may influence 
the properties of the product surface. ** Thus, it is critical to under¬ 
stand the rates of adsorption/desorption phenomena taking place at 
the interface between a solid and a polymeric melt. 

In addition to the study of polymer melt dynamics near flat walls, 
there is a growing interest in the direction of studying nanopar¬ 
ticles dispersed in a polymeric matrix. This is ultimately re¬ 
lated to the industrial importance of nanocomposites, in which the 
filler particles may have dimensions comparable to those of poly¬ 
meric chains. *^ Several simulation studies have appeared recently 
in which static and dynamic properties of model nanocompos¬ 
ites were investigated. *^“** Alongside with molecular simulations, 
mesoscopic models invoking the picture of “glassy bridge” forma¬ 
tion around the dispersed nanoparticles have been developed. *®’^® 
Models of this kind assume that the polymer dynamics around the 
dispersed particles are orders of magnitude slower that in the melt. 
The “percolation of free volume distribution” model^* was used 
by Berriot and co-workers for explaining the microscopic ori¬ 
gin of the reinforcement in filled elastomers, as a consequence of 
the presence of a gradient in glass transition temperature around 
the fillers. Estimates of polymer adsorption/desorption rates at the 
solid surface may be of help in parametrizing these models in a 
bottom-up approach. 

Erantz and Granick^'* were the first to study experimentally the 
adsorption-desorption kinetics of polymer chains from solution to 
a solid surface. They monitored, by in situ Fourier-transform in¬ 
frared spectroscopy, the surface coverage of polystyrene adsorb¬ 
ing to oxidized silicon from cyclohexane. These authors have dis¬ 
tinguished three time scales: for diffusion to and adsorption onto 
an initially bare surface (which is a rapid procedure), for ensuing 
surface rearrangements (slow and history dependent), and for dis¬ 
placement. The duration of the last part (i.e. the displacement of 
adsorbed chains along the solid surface) grows exponentially with 
the molecular weight of polymeric chains. Even at overall equilib¬ 
rium there is a steady-state traffic of macromolecules between the 
adsorbed and the unadsorbed state. 

Despite the fact that plenty of simulations have been conducted 
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on polymer/solid systems, only few of them have quantified ad¬ 
sorption/desorption kinetics. Smith et al.^^ have carried out molec¬ 
ular dynamics (MD) simulations of a coarse-grained bead-spring 
model in order to examine the static and dynamic properties of 
polymer chains in a melt in the presence of a solid surface. These 
authors have shown that the population of chains contains a sub¬ 
set of weakly adsorbed chains with rapid desorption dynamics and 
strongly adsorbed chains with slower desorption dynamics. The 
strongly adsorbed chains were held responsible for a slow-down of 
dynamics near the surface. The desorption time was found to scale 
slightly more slowly than (exponent 1.5 - 2) with N being the 
length of the chain. The observation holds even for chains that are 
in the reptation regime (where scaling could be expected). If 
only simple diffusion of the adsorbed beads were taking place, a 
kinetic model employed in that study would yield a desorption time 
scaling as N'^. The MD results indicate a more complex behavior 
of chains that does not depend only on the number of adsorbed seg¬ 
ments, but probably also on the time span for which a chain has 
been attached to the surface (which is inversely related to the dis¬ 
tance of the chain center of mass from the surface). 

Yelash et al.^^ have employed MD simulations of a chemi¬ 
cally realistic model of 1,4-polybutadiene between graphite walls. 
They have shown that attractive confining walls introduce one more 
mechanism (in addition to the caging effect and intramolecular 
conformational barriers) for time scale separation of glass form¬ 
ing polymers. This mechanism is the slow desorption kinetics of 
monomers from the surface, leading to slow layer exchange dy¬ 
namics on the length scale of (at least) the radius of gyration of the 
chains. A layer-wise analysis of the relaxation times, as a function 
of the distance from the solid surface, has revealed that the adsorp¬ 
tion/desorption kinetics of monomers at the wall,^^’^® give rise to a 
two-step relaxation behavior, as the other two mechanisms also do. 
In a later study, they have shown that a three-step (one for every 
timescale separation mechanism) decay can be observed in inco¬ 
herent scattering experiments and discussed its relevance for the 
glass transition of confined polymers by analogy to surface critical 
phenomena. 

One of the most common approaches used for the study of equi¬ 
librium properties of polymer/solid interfacial systems is the self- 
consistent field (SCF) method. The essence of the SCF method 
is the replacement of the ensemble of interacting polymer chains 
with a system of noninteracting chains subject to some position- 
dependent (complex) chemical potential fields. These fields dictate 
the conformations of the polymer molecules, imposing the spatial 
distribution of the polymer. On the other side, the chemical poten¬ 
tial fields depend on the polymer density distribution. Thus, the tar¬ 
get of this approach is to determine the fields in such a way that they 
are consistent with the spatial distribution of polymer they create. 
A very important step in solving the SCF problem is the identifica¬ 
tion of stationary field configurations that correspond to extrema of 
a complex effective Hamiltonian. For a real field theory, such a 
configuration can be a local minimum, maximum, or a saddle point 
in the field configuration space. However, when the field variables 
are extended to the complex plane, the energy surface is generally 
saddle shaped in the vicinity of an extremum, so we shall refer to 
stationary field configurations as saddle points. It has been proved 
that there is an explicit connection between the saddle point chem¬ 
ical potential field configuration of the SCF problem and a corre¬ 
sponding mean-field approximation to the free energy. 

There are two general classes of techniques that have been ap¬ 
plied to solving SCF equations and hence to the numerical com¬ 
putation of saddle points. The first are spectral methods, which 
attempt to represent the various spatially varying fields in a trun¬ 
cated Fourier-like basis. Alternatively, the nonlinear equations can 


be tackled in real space by suitable finite difference or finite element 
discretization of a computational domain. Daoulas et al.^^ have 
determined the equilibrium properties of a polymer melt of specific 
chemical constitution (polyethylene), adsorbed on a certain solid 
substrate (graphite), through a continuum SCF approach. While 
the monomer-monomer interactions were introduced in the simpli¬ 
fied Helfand approximation, interaction potential between 

the polymer and the substrate was tuned in order to reproduce the 
total energy of adsorption and the characteristic length scale of the 
density variations of the atomistically simulated system.^* In con¬ 
trast to previous SCF studies, for the first time, the chain connec¬ 
tivity was represented through the wormlike chain model, in addi¬ 
tion to the simpler Gaussian model. Polymeric chains were found 
significantly “flattened” near the boundary, while perturbations in 
their conformational properties due to the surface persisted over 
a characteristic length which is roughly 1.5 times the unperturbed 
mean radius of gyration. The conformations of adsorbed molecules 
were further characterized by considering the properties of tails, 
loops and trains, as derived by the SCF analysis of the two mod¬ 
els (wormlike and Gaussian) and the atomistic simulations. SCF 
predictions from the wormlike chain model were found closer to 
the atomistic simulation results for short chain systems. The supe¬ 
rior performance of the semiflexible chain model when compared to 
the Gaussian one, is due to the introduction of an additional length 
scale, that of the persistence length. However, as the molecular 
length increases, the effect of local chain structure on global chain 
conformational behavior decreases, rendering the differences be¬ 
tween the two models insignificant. 

In a recent study, Klushin et al."*^ have investigated the effects of 
the range of adsorption potential on the equilibrium behavior of a 
single polymer chain end-attached to a solid surface. The exact an¬ 
alytical theory for ideal lattice chains interacting with a planar sur¬ 
face via a square-well potential was presented and compared to con¬ 
tinuum model results and to Monte Carlo (MC) simulations using 
the pruned-enriched Rosenbluth method for self-avoiding chains on 
a simple cubic lattice. We are particularly interested in the analyt¬ 
ical approach, where the partition function of a chain with both its 
ends fixed in space is obtained as the solution of the partial differen¬ 
tial equation of diffusion of the Green’s function, an argument we 
will also invoke in our methodology. In that framework, the adsorp¬ 
tion order parameter (the bound fraction of monomers) was found 
to be a function of a single argument, the potential well depth. Tail, 
loop, and train distributions at the critical point were evaluated by 
MC simulations and compared to analytical results for ideal chains 
and with scaling theory predictions. The behavior of a self-avoiding 
chain was found to be remarkably close to that of an ideal chain in 
several aspects. 

The present work is aimed to address the kinetics of adsorption 
and desorption of an end-constrained polymer strand belonging to 
a melt close to an adsorbing solid substrate. In real-life polymer 
systems, end constraints may be created by crosslinks or entangle¬ 
ments with other chains and may fluctuate thermally. Such thermal 
fluctuations of the end positions are not considered here, in order to 
facilitate a clean mathematical formulation. The ends of the con¬ 
strained chain, which in practical situations may be a subchain of 
a longer chain or of a three-dimensional network, are considered 
as fixed in space. Fluctuations in the chain ends can be introduced 
through a straightforward generalization of the formalism presented 
herein. We employ the “Gaussian thread model”, which considers 
the polymeric chains as fully flexible and infinitely extensible, con¬ 
tinuous threads or “paths” and offers the advantage of facilitating 
the analytical or numerical evaluation of the chain conformations 
through the solution of the Edwards equation in the mean field ap¬ 
proximation. The fields obtained from the solution of the SCF 
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problem are then used in order to calculate the partition function 
of end-constrained chains, as a function of the distance of their 
ends from the solid surface. This allows us to define a potential of 
mean force experienced by an end-constrained chain as a function 
of its distance of closest approach to the solid surface. The fea¬ 
tures of the one-dimensional free energy profile are incorporated in 
a dynamically-corrected transition state theory (TST) formalism in 
order to estimate adsorption and desorption rate constants, as func¬ 
tions of the molecular weight of the chain and the distances of its 
ends from the substrate. These rates have been validated against 
rates obtained by hazard-plot analysis of MD trajectories of the 
same system."*^ The knowledge obtained from this study will be 
of great value for developing mesoscopic simulation approaches to 
treat rate-dependent deformation and flow near solid walls or in the 
presence of nanoparticles. We expect that the characteristic time as¬ 
sociated with exchange dynamics (adsorption/desorption) of chain 
strands from the melt to the surface may play an important role in 
the enhancement of the mechanical properties of filled elastomers. 

The system considered in this work is a polyethylene (PE) melt 
confined between two graphite phases. PE/graphite nanocom¬ 
posites are of paramount industrial importance because they 
exhibit electrical and thermal conductivity,^® holding a great 
promise as antistatic materials, capable of dissipating static charges 
safely from component surfaces. The thermal conductivity of 
PE/graphite composites can be as high as 11.28 
These composites exhibit a temperature-dependent electrical con¬ 
ductivity due to thermal expansion-dependent percolation of the 
carbon particles.The electrical resistivity of carbon black 
filled with high density polyethylene increases significantly when 
the composite is heated to the melting temperature of the matrix. 
However, void formation was observed to occur during polyethy¬ 
lene crystallization as a result of filler particle/polymer matrix in¬ 
teractions, by small angle neutron scattering. Moreover, PE is not 
too different from cis-1,4 polyisoprene and the surface of carbon 
black could be approximated as graphite, so our system is close, at 
a molecular level, to carbon black-reinforced rubbers. 


SCF Formulation for a Homopolymer Melt Next to 
an Adsorbing Solid Surface 

Field Theoretic Formulation of the Grand Partition 
Function 

The continuum SCF approach is based on a path-integral represen¬ 
tation of the partition function. In the past, this approach has been 
extensively employed in numerous works. Here, an effort 

will be made to present the theory as briefly as possible, empha¬ 
sizing only essential points specific to the present work. We con¬ 
sider a polymer melt in the presence of a solid phase, which exerts 
a position-dependent field Us (r) on the polymer segments. Each 
chain is envisioned as a sequence of N statistical segments obeying 
Gaussian statistics, and Us (r) is an energy per statistical segment 
located at r. We focus on a region of total volume V occupied 
by polymer at temperature T. The boundaries of that region are 
partly defined by the surface(s) of the solid. The polymer in the 
considered control volume is at equilibrium with a bulk polymer 
phase, with which it can freely exchange chains. The chemical po¬ 
tential of a chain in that phase is flN (where by definition, fl is 
the chemical potential per segment). The polymer in the interfacial 
region conforms to the probability distribution of the grand canon¬ 
ical ensemble. Following the statistical mechanical development 
of Daoulas et al.,^^ based on a simplified, Helfand-type effective 
Hamiltonian for polymer-polymer interactions, the grand canonical 


partition function can be written as in Eq.(2) of that paper: 


S =£ 1 exp f / n (•) ^ [>•« (•)] 

-^y't/s(r)0(r)dV) (1) 

where is the Boltzmann constant, po and Kj are the mean seg¬ 
ment density and isothermal compressibility in the bulk melt, while 
p, N, and n are the segment chemical potential, the chain length 
of the coarse-grained Gaussian thread, and the number of chains, 
respectively. M is a normalizing prefactor, and (•) denotes 
functional integration over all possible conformations (“paths”) of 
chain a, while 0 (r) is the volume fraction operator, given by 


/ ^(r-ra(.s))d^ 

0 (r) = - ( 2 ) 

Po 

where 5 is a scaled variable measuring how far along the contour 
length of the chain a considered segment lies; it ranges from 0 
(chain start) to 1 (chain end). The functional ^ [ra (•)] accounts 
for the chain connectivity, and in the case of the Gaussian model it 
takes the form 


^ [ra (•)] = exp 



(3) 


where is the mean-squared radius of gyration of the polymer 
chain in the bulk. It can be seen that, in the case of the Gaus¬ 
sian model, only one characteristic conformational parameter of the 
atomistic polymer chain enters the coarse grained model: R^. 

Following a standard field theoretical approach, it is possible to 
replace the system of interacting chains subject to the field of the 
solid with a system of noninteracting chains subject to the field of 
the solid and an additional fluctuating field representing the remain¬ 
ing chains. This procedure is described in detail in Appendix A. We 
introduce the notation Q [iw + t/s] to indicate the partition function 
of a single chain subject to the field iw -|- Us acting on its segments, 
relative to the partition function of a field-free chain. The grand 
partition function can be written as: 

E = cJ &[j5w]exp{—l5H[w]} (4) 

with 

H M =J [-ipow (r) -f y (pow 

- J exp ^ j ^ZfreeG K + Us] (5) 


It is remarkable that the the effective “Hamiltonian” H [w] incorpo¬ 
rates a term proportional to exp ^ ^ and to Q, rather than to pA 

and to In 2, because of the grand canonical formulation adopted. 
Zfree = / (•) [fa (•)] is the partition function of a free chain. 

By invoking a saddle point approximation (which is described in 
detail in Appendix B), the effective Hamiltonian, Eq. (5), at the 
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optimum can be written as: 


H = -poJ dVw'(r) 0 (r) + dV[l - (r)]^ 

+Po J dV[/s (r) 0 (r) 

-^exp (|^)AfZfree2[iw + i7sl (6) 

where we have set w' (r) = i w (r) + Us (r), a real field (see Ap¬ 
pendix A). The effective field w' (r) and the volume fraction, (j) (r), 
are coupled to each other through a set of equations depending on 
the model adopted for the chain connectivity. For the case of the 
Gaussian threads considered in this work, assuming that distances 
between solid surfaces present are large in comparison to Rg, such 
that bulk conditions prevail far from the surfaces in the interfacial 
polymer, these equations are formulated as: 

w' (r) = [0 (r) ~l] + Us (r) (7) 

KtPo 

which is the self-consistent field expression for the real field and 

^{r)=f dsq{r,s)q{r,l-s) (8) 

Jo 

which stands for the segment balance condition. In the formulation 
adopted here, (j) (r) is a segment density reduced by its value po in 
the bulk melt at the same polymer chemical potential. The quantity 
q (r, i) is the restricted partition function, proportional to the prob¬ 
ability density that the segment at fractional contour length i- of a 
chain which may have started anywhere in the system lies at po¬ 
sition r. It is normalized by the corresponding probability density 
in the bulk melt and obeys the initial condition ^(r,0) = 1. The 
restricted partition function ^(r,i') can be calculated through the 
Edwards diffusion equation: 

='^gV??(r,i)-j8A>v'(r)^(r,j) (9) 

A relation exists between the restricted partition function q{r,s) 
and the single chain partition function in the presence of the field, 
Q[iw + Us] = QW]-. 

Q[w']=^ j d\q{r,l) (10) 

The SCF problem to be solved in order to capture structure 
and thermodynamics in the interfacial system consists of the par¬ 
tial differential “diffusion” equation, Eq. (9), the definition of the 
the self-consistent field, Eq. (7), and the segment balance Eq. (8). 
These must be solved numerically in the unknown restricted parti¬ 
tion function q (r, j), the volume fraction profile (j) (r), and the self- 
consistent field w' (r). As already pointed out, the “initial” condi¬ 
tion ^ (r, 0) = 1 applies. Additionally, boundary conditions 

must be used that reflect the geometry of the problem. 


main 0 < z < and 0 < .s < 1: 
dq{z,s) _ 


(/i(z) = / dj^(z,i)^(z,l-j) 

Jo 

W(z) = ^W{z)-l]+^NUs{z) 
KtPo 


under the boundary conditions: 

dq(z,s) 

dz 

and the initial condition: 


q{o,s) = o 


= 0 




( 11 ) 

( 12 ) 

(13) 

(14) 

(15) 


^(z,0) = 1, Vz > 0 (16) 

The plane of symmetry at z = Eq. (15), implies that the poly¬ 
mer is sandwiched between two identical surfaces. The problem 
is solved within half of the gap only. We are here interested in 
cases where the thickness of the interfacial domain, in which the 
problem is solved, is large enough, in comparison to the range of 
the segment-surface interactions and in comparison to the chain di¬ 
mensions Rg, for bulk conditions to prevail at z = L^- We expect 
that Us (iz) = 0 and (j) (L^) = 1, hence W (L^) = 0 as well. How¬ 
ever, it is certainly possible to address thin films of polymer, such 
that bulk conditions are not realized anywhere in the film. 

Commonly, 57-59 considering questions related to polymer 

adsorption, the Edwards diffusion equation is supplied with an ef¬ 
fective boundary condition. This practice originates from the work 
of de Gennes and amounts to splitting the total effective field felt 
by a polymer segment into a part originating from its interaction 
with the substrate and a mean-field part accounting for the interac¬ 
tions with the rest of the polymer segments. In our case, however, 
following Daoulas et al.,57 we have to preserve the surface com¬ 
ponent of the effective field, W (z). To this end, the propagator 
equation, Eq. (11), has to be supplied with an absorbing boundary 
condition, Eq. (14), requiring q{z = 0,i) = 0. Wu et al.®* gave a 
detailed reasoning of how, beyond some distance from an enthalpi- 
cally neutral substrate, the microscopic absorbing boundary condi¬ 
tion effectively reduces to the well-known for this case reflecting 
boundary condition. 

The SCF equations formulated in the current work (i.e. eqs 
Eq. (11), Eq. (12) and Eq. (13)) can be solved following a Newton- 
Raphson method, which is expected to converge fast, albeit, requir¬ 
ing more memory and CPU time per iteration than the more com¬ 
monly used successive substitution, due to the need to compute and 
store the Jacobian. The details of the discretization and numerical 
solution of the system of eqs Eq. (11) -Eq. (16) can be found in the 
Supporting Information to the present paper. 


Numerical Solution of the SCF Problem 

We specialize to the one-dimensional problem, where the poly¬ 
meric melt is conflned between two identical planar surfaces, within 
a gap of width 2L^. After recalling the translational invariance of 
the system along the xy plane, we can write q(r,s) = q{z,s) and 
(j) (r) = (j) (z). We have to solve the following integrodifferential 
system in the unknown functions q (z,i), (/> (z) and W (z) in the do- 


SCF-Based Formulation for Estimating the Rates 
of Adsorption and Desorption of Entangled or 
Crosslinked Strands from a Melt or Rubbery Poly¬ 
mer 

Eormulation of Adsorption and Desorption as Infre¬ 
quent Event Processes 

We consider a polymer above its glass transition temperature, con¬ 
sisting of randomly coiled strands whose ends are more or less con¬ 
strained by entanglements or crosslinks, next to a flat solid surface 
exerting a potential energy field Us (r) on polymer segments. Our 
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objective is to quantify the rates of adsorption and desorption of a 
strand in the sea of other strands, as functions of the positions of 
the two ends r(^ = 0) and r (i = 1), with s being the reduced con¬ 
tour length along the strand. In the following we will use the terms 
“strand” and “chain” interchangeably, even though the considered 
strand may in fact be a subchain of a longer chain. Without loss of 
generality, we consider within a linear monodisperse melt a chain 
of length N statistical segments whose ends are constrained at or 
near r(0) and r(l). The conformations of such a chain in the in¬ 
terfacial region can be analyzed via a SCF model of the type we 
discussed above. This model provides probability distributions for 
the ends (i- = 0 and i = 1) to find themselves in specific positions 
in space. 


Z 

-z = 0 


Figure 1. Schematic representation of an adsorbed chain. In this example, 
one of the chain ends, r(0), is a crosslink, while the other one, r(l), is an 
entanglement point with another chain (marked in black). The height of the 
layer, in which adsorption is considered to take place, is S. 

Based on the form of the segment-surface interactions, Us{z), 
one can define a thickness 8 next to the surface such that, if a seg¬ 
ment lies at 0 < z < 5 with z = 0 indicating the plane of the surface, 
it is considered as adsorbed, while if it lies at z > 5, it is considered 
as free. A chain will be considered as adsorbed (subscript “a”) if it 
contains at least one adsorbed segment. Otherwise it will be con¬ 
sidered as free (subscript “f”). An adsorbed chain, lying between 
a crosslink and an entanglement, is depicted in Figure Figure 1, 
along with the basic geometric quantities used in the formulation. 
The chain may execute a very complex motion near the surface, in¬ 
volving fluctuations in the number and length of trains, loops and 
tails.In this work, we intend to derive overall rate constants for 
the rates of adsorption and desorption of the end-constrained chain 
on the surface, which are applicable to the long-time behavior of 
the chain and are integrated over all these fluctuations. With the 
chain ends fixed at r(0) and r(l) in the z> 8 domain, we think of 
adsorption as a transition between a “free” state (z(j) > 8) and an 
“adsorbed” state (z(x) < 8 for some j in 0 < x < 1). 



Configurational Integral of End-constrained 
Chains 


Thus, the propagator G>j (r,-, r', x) becomes: 

G>z (r;,r',i) = Gx Gy {yi,y',s) G>z {zi,z,s) (17) 

with the factors Gx(x/,x,j), Gy{yi,y,s) being unconstrained one¬ 
dimensional Gaussian functions emanating from x;, y;, respectively, 
and therefore independent of z. By definition. 


G>z [zi,z!,s) = 0 Vz' forz > z; (18) 

so G>j (z^,z^^) should be computed for 0 < z < Zi and ^ > z. 

The configurational integral Z(z) is related to G>j (z;, z^ j) as fol¬ 
lows: 


^bulk 


= [ dz' r dx r dy f 

Jz'>z 7 — 00 7 — 00 Jo 


ds 


[G>z (z(0) , z' , i) G>j (z( 1), z' , 1 - .s) 
X Gx (x(0),x,x) Gx (x(l),x, 1 - s) 

X Gy (y (0), y, x) Gj, (y(1), y, 1 - j) ] 


(19) 


where .^ulk is the configurational partition function of a free chain 
of length N starting at a given point in the melt. It should be 
noted that 3j,uik is independent of the exact starting point of the 
free chain. Note also that the ratio 2f{z)/^hu\k has dimensions of 
inverse volume, since i2^ulk in the denominator has only one end 
constrained, while ^ (z) has two ends constrained. The right-hand 
side of Eq. (19) also has dimensions of inverse volume, since the 
propagators along all x, y and z have dimensions of inverse length. 
The propagator G>J■(Zi,z^x) needed within Eq. (19) for z; = z(0) 
and Zi = z(l) satisfies the following differential equation: 

<5G>j(zi,z',i) _ 2^^G>j(z/,z',i) 
ds S ^ 2/2 

-W(z')G>,(zi,z',^),i>0 
G>z {zi,z!,Q) = 8 (z! —Zi) (initial condition) 

G>j (zi,z,0) = 0 (boundary condition at z^ = z ) 
lim G>z izi,z',s) = 0 (boundary condition at z' °°) (20) 

7 ' —ioo 


The absorbing boundary condition at z^ = z stems from the require¬ 
ment that the entire contour of the considered chain must lie at 
z' > z. The self-consistent field VF (z') is already known from our 
solution to the SCE problem. Eq Eq. (20) must be solved for given 
Zi, z with 0 < z < z; to determine G>j {zi,z!,s) for all z' > z. 

Eor our one-dimensional problem, the integrations over x and y 
appearing in Eq. (19) can be carried out analytically, yielding: 


I dx'Gx(x(0),x',x)Gx(x(l),x',l-y) 
1 [ [x(l)-x(0)f 

= 2^“P — 47^2 - 


( 21 ) 


In the context of the continuous SCF model of a melt next to an ad¬ 
sorbing solid surface, we define 3f{z) as a configurational integral 
over all chains, with their ends fixed at r(0) and r(l), which do not 
have any segments in the region 0 < z' < z, with z' = 0 denoting 
the substrate surface and z' = 8 denoting the dividing surface for 
adsorption. (z) is expected to be a decreasing function of z. 

For the calculation of 2f{z) we define a propagator G>j (r,',r',j) 
for a chain that has started at position r, to be at position r' at con¬ 
tour length s, the whole contour of the chain from 0 to s lying en¬ 
tirely in the region of space z' > z. This is equivalent to posing an 
imaginary neutral solid wall at z, while solving the diffusion dif¬ 
ferential equation, Eq. (20). We expect G>j. (r,',r', j) to be factoriz- 
able, due to the infinite extensibility of the Gaussian chain model. 
Presumably, for finite-contour length models, like the freely-jointed 
or the worm-like chain models, such a factorization is not possible. 


independent of s, as one would expect, since we have an uncon¬ 
strained Gaussian chain in the x-direction. Eq Eq. (21) is readily 
applicable in the y-direction, too. Erom eqs Eq. (19) and Eq. (21): 

nz) \ [x(l)-x(0)p + [y(l)-y(0)f ' 

^bulk P [ 4Rl \ 

[ dz' f'dsG>z{z{0),z',s)G>z{zW,z,l-s) (22) 
Jz'>z JO 

Let and be the configurational integrals of the chain in the 
adsorbed and free states, for given r(0) and r(l). The former one 
can be cast as the difference between the configurational integral at 
z = 0 (all chains having their segments closest to the surface above 
the z = 0 level) and the configurational integral at z = 5 (all chains 
having their segments closest to the surface above the z = 8 level). 
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With these definitions: 


a Boltzmann distribution): 


jrf = jr(S),jr^ = ^(o)-jr(S) (23) 


and the equilibrium probabilities Pf and P-^ of the chain being free 
and adsorbed, respectively, are: 


_ iTf „ _ 

iTf + J^a’ ^ iTf + ^a 


(24) 


Connection of Infrequent Event Formulation with 
SCF Model 


Clearly, associated with the configurational integral S’ (z) is a po¬ 
tential of mean force (configurational free energy) for a chain to 
have its segment(s) lying closest to the surface at z. We envision 
that, in order to desorb from the adsorbed state or adsorb from the 
free state, the chain must pass through a “transition state” configu¬ 
ration that has at least one segment in the plane z = 5 but all other 
segments at z > 5. In other words, we consider as “reaction coor¬ 
dinate” the lowest among all z-coordinates of the segments of the 
considered chain and a “dividing surface” for this reaction coordi¬ 
nate at z = 5. By definition, then, —diF(z)/dz is a configurational 
integral per unit length along the z-direction for a chain, with its 
ends fixed at r(0) and r(l), to have its segment(s) closest to the 
surface between z and z + dz. Thus, the configurational integral per 
unit length at the dividing surface is: 

= - — (25) 

while the potential of mean force can be obtained as: 


A(z) = 


dS’(z) 

dz 


C4 


(26) 


where a multiplicative constant with the units of length to the fourth 
power, C 4 , is introduced in the potential of mean force, in order to 
make the argument of the logarithm dimensionless. The introduc¬ 
tion of the constant shifts the reference level for the estimation of 
free energies without affecting their physical significance. We ex¬ 
pect the potential of mean force, A(z), to exhibit a local maximum 
with respect to z, which will be used as an estimate of the position 
of the dividing surface. 

The next step in calculating the rates of adsorption and des¬ 
orption would now be to apply dynamically corrected transition 
state theory, based on the theory of Bennet® and Chandler.®^ It 
would be reasonable to assume that the chain is well described by 
Kramers®^ theory in the high friction limit. If we consider a Rouse 
picture for the chain in the melt, and let ^strand denote a friction co¬ 
efficient between a strand and the rest of the melt (proportionality 
constant between frictional force exerted by the melt and velocity 
of the strand relative to the melt, measured in units of mass per 
time), the dynamical correction factor associated with recrossing 
the dividing surface will be: 


/ = 


\/ ^b^strand 


Cstrand 


(27) 


where Kj, is related to the curvature of the free energy profile at the 
top of the barrier: 


K'b = - 


d^A(z) 

dz^ 


z=S 


(28) 


By assuming an effective mass mjtrand for a polymer strand, at 
temperature T the strand will have a mean thermal speed along z 
(mean of the magnitude of the velocity in one dimension, assuming 


(kzl) 


2kYiT 
^^strand / 


(29) 


Thus, the adsorption rate constant for the end-constrained chain 
emerges as: 


1 

^f-^a ~ 2 ^ (kzl) 


1 1 


1/2 j 

1 1 

f dS 


Cstrand 

V 2 ;r ) 

' S{5) ' 

V dz 


(30) 


with iq, being given in terms of S’ (z) via eqs Eq. (28) and Eq. (26). 
Similarly, the desorption rate constant for the end-constrained chain 
emerges as: 


1 

^a^f (kzl) 

1 1 

fstrandV 27r ) S (Q) - S (8) 


dS \ 
dz z=s) 


(31) 


Connection to Atomistic Simulations 


Representation of the PE Melt/Graphite Interfacial 
System in the SCF/TST framework 

After the set of SCF and TST equations has been formulated, there 
are several parameters to be determined, so that a correspondence 
between the field-theoretic and the atomistic PE/graphite systems 
can be established. Since all atomistic simulations were performed 
at T = 450 K, the same temperature will be assumed during the 
SCF/TST calculations of the present work. The bulk density of 
the polymer is po = 0.766 g/cm^. The compressibility appearing 
in the SCF equations, (fj, will equal the isothermal compressibil¬ 
ity of the atomistically studied bulk polymer.Taking into ac¬ 
count our atomistic simulations and experimental data,®^ we em¬ 
ploy iCp = 1.43 X 10^® Pa^ *. The length scales of the field-theoretic 
and the atomistic molecular models are related through two invari¬ 
ant parameters: the contour length of the chain, L, and the mean 
squared radius of gyration of the chain, R^. The latter can be calcu¬ 
lated as Rg = {C/6)N\,l^ where C is the polyethylene characteristic 
ratio, A(, the number of chemical bonds along the chain and /b their 
bond length. The bond length used in the atomistic simulations 
is /b = 1 -54 A and we set C for the molecular weight of interest, 
following the work of Karayiannis et al. A crucial step in imple¬ 
menting the field-theoretic representation for the PE/graphite sys¬ 
tem is the definition of the polymer/substrate interaction potential, 
(/s(z)- In this scope, we follow the work of Daoulas et al.^^ where 
a square-well potential was used. The depth of the square well po¬ 
tential was determined by requiring that the field-theoretic model 
have the same adsorption energy per unit surface as the atomistic 
one. 

The Rouse model considers a polymer chain as a sequence 
of Brownian particles connected by harmonic springs and mov¬ 
ing in a viscous medium representing the background environment 
formed by all other chains. An important parameter in the Rouse 
formulation is the friction factor §3 (measured in units of mass 
per time), characterizing the resistance encountered by a monomer 
unit moving through its surroundings. It can be thought of as 
the proportionality constant between the velocity of the Brown¬ 
ian particle and the frictional force exerted on the particle as it 
moves through the “sea” formed by all other particles in the sys¬ 
tem. Harmandaris et al.®* have estimated the segmental friction 
coefficient of polyethylene carbon atoms by mapping atomisti- 
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cally detailed MD simulations of linear PE melts onto the Rouse 
model. Their simulations have shown that there is a minimum 
chain length (around Cgg), above which i^o can be considered as 
constant, chain length independent parameter of the system. The 
asymptotic value is ~ 4.5 x 10^*^ kg/s. If we think of the fric¬ 
tion coefficient as being proportional to the mass of the entity it 
refers to, we can estimate the maximum possible friction coeffi¬ 
cient for the diffusion of the strands treated in the TST framework 
as: ^strand = (™strand/'ncH 2 ) Co- The real friction coefficient of the 
strand will lie between the extreme values: ^3 < ^strand < Co m'c"*' ' 
We do not take into account the possible dependence of the fric¬ 
tion coefficient fg on the distance of the strand from the graphite 
surface. 

Atomistic Simulations 

The atomistic simulations were conducted by describing the PE 
chains with a fully flexible united atom model, which considers 
each methylene (CHj) and methyl (CHj) group along the chain 
backbone as one interaction site.^*’69 polymer intramolecular in¬ 
teractions (bond stretching, bond angle bending, and torsional po¬ 
tential) are implemented following Nath et al. Non-bonded inter¬ 
actions between pairs of CH 2 or CH 3 are described by a pairwise- 
additive Lennard-Jones 12-6 potential with input parameters pro¬ 
vided by the TraPPE forcefield of Martin and Siepmann. Ewald 
summation has been applied to the attractive part of all polymer- 
polymer nonbonded interactions, in order to avoid distortions 
of the density profile and associated stress profile. To calculate 
the potential energy of interaction between a polymer atom and a 
semiinfinite graphite phase, we implement an efficient and accu¬ 
rate Eourier summation method designed by Steele, which takes 
advantage of the symmetry in the crystalline substrate. Eor more 
details and a thorough study by using atomistic simulations the in¬ 
terested reader is referred to our accompanying work on the same 
system. In the following paragraphs we will summarize the vital 
pieces of information relative to the present work. 

Initial configurations of polyethylene confined between two 
graphite phases have been obtained by a MC builder capable of 
building polymeric chains of arbitrary geometry in heavily con¬ 
strained environments, based on the quasi-Metropolis scheme of 
Theodorou and Suter. Then, connectivity altering moves, such 
as concerted rotationand double bridging®® were employed 
in MC equilibration of the initial configurations. The simulations 
were carried out in orthorhombic cells with periodic boundary con¬ 
ditions along the x and y directions, while the system along the z 
direction was considered as finite. All MD simulations have been 
conducted using the Large-scale Atomic/Molecular Massively Par¬ 
allel Simulator (LAMMPS),^^ extended with the graphite potential 
of Steele, which we have incorporated into the LAMMPS source 
code. 

In our MD simulations, the graphite substrates lie at the z = 0 
and z = TLj planes (Figure Figure 2). The simulated polymeric 
melt consists of 79 chains of 200 united atoms (either CH 2 or CH 3 ) 
each. All simulations (MC and MD) have been carried out at the 
temperature of 450 K and lateral pressure = 1 atm. After letting 
the system relax in the course of an MD simulation, we decide to 
anchor an internal strand (marked in red in Figure Figure 2) by ze¬ 
roing the forces acting on its ends, having ensured that the velocities 
of its ends are also zero. We then monitor the residence time of this 
strand in the free (all united atoms at z > 5) and adsorbed (at least 
one united atom at z < 5) states. The length of the runs was in the 
order of 400 to 500 ns, by employing the RESPA multiple-timestep 
integrator^* with the outer integration step being 4 fs. 

The adhesion tension, i.e., the difference between the surface 



Figure 2. Schematic illustration of the end-constrained subchain of 50 
CH 2 units (marked in red), whose residence times in the adsorbed and des¬ 
orbed state, in the course of an MD simulation, were used for estimating 
adsorption and desorption kinetics via hazard-plot analysis. The chain is 
depicted as a necklace of the van der Waals spheres of the atoms, with atoms 
not belonging to the selected subchain shown in blue. The ends of the red 
subchain have been anchored at height z(0) = z(l) = 1.3 nm. The remain¬ 
ing chains in the system are depicted in wireframe representation, for clarity. 
Two indicative layers of graphite atoms lie at the upper and lower z bounds 
of the simulation box. Visualization made by using the VMD software. 


tension of the clean graphite surface, 7 s, and the interfacial ten¬ 
sion of the graphite/polymer interface, 7 fs, can be computed follow¬ 
ing the approach developed by Tolman and refined by Kirkwood 
and Buff, as an integral of the ensemble averaged instantaneous 
difference between the normal and tangential stresses CTx(z) and 

(Ts - Tfs) = - ^ ' (tT|| (z) - (z)) dz^ (32) 

where, in our geometry, CTx(z) = dzziz) and 0 '||(z) = 
(cTj:;c(z) + CTy.y(z))/2.*^ Awuy from an interface, CT^ = and 
the integrand vanishes. Therefore the nonzero contributions to the 
integral in Eq. (32) come from the interfacial region. The outer fac¬ 
tor of 1 /2 in Eq. (32) accounts for the presence of two liquid-solid 
interfaces. The stress tensor as a function of the distance from the 
surface can be estimated by summing atomic-level stress contri¬ 
butions, CTj, of atoms lying in a specific slab along the z direction, 
^ (z) = (1/Tslab (z)) LiGslab Here, we use a Voronoi tessel¬ 

lation to define atomic volumes,®®’^® such that summing over all 
atoms gives the total volume of the system, V = £,4/. Following 
Mansfield and Theodorou and Theodorou et al. we define the 
atomic level stress as: 


C/.LM = 


1 


1 

2Vi 


J^i 


\ min.im. 

Oi) 


irmin.im. 

fij,M 


1 WwaJls 


(33) 


where r,- and are the position vectors of atoms i and j, v; the 
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velocity of atom i and F/y is the force exerted on atom i by atom 
j. In general, the force F;j between two sites i and j is defined as 
F,y = ~Vr,-r .l^ where Y is the total potential energy of the sys¬ 
tem and the gradient is taken keeping all intersite separation vec¬ 
tors other than r, — constant. The superscript “min.im.” indi¬ 
cates interatomic distances and forces calculated according to the 
“minimum image convention”. The indices L and M, indicating 
the three coordinate directions in a Cartesian system, assume the 
values X, y and z. The last term of Eq. (33) stems from the interac¬ 
tion of the atoms with the substrate. The quantity dik,L multiplying 
the force exerted by the k-th wall, on atom i, takes the val¬ 
ues: = Xi, dik^y = yt and dik^z = “ z*:,w with Zk^vj indicating 

the position of the k-th (upper or lower) graphite slab along the z- 
direction. Moreover, long-range contributions to the atomic-level 
stresses of Eq. (33) have been calculated based on eqs 27 and 28 of 
ref 72. 

Results 

Self-Consistent Field and Local Volume Fraction 

The SC field obtained from the solution of the Edwards diffusion 
equation ( Eq. (9)) is presented in Eigure Eigure 3, for two different 
chain lengths, N = 50 and N = 100. Along with the field, the attrac¬ 
tive potential of the solid surface is also shown for comparison. The 
SCE fields exhibit a smooth behavior due to the fine discretization 
used for the solution of the SCE problem. The square well form 
of the potential causes a minor numerical instability at the distance 
where it becomes zero, yielding an discontinuous behavior of the 
SCE profiles,at the end of the well. This is expected and does not 
affect the following calculations. The use of a square-well potential 
is in the spirit of coarse-graining, given the fact that this potential 
preserves the same adsorption energy per unit area, and the range of 
the polymer density fluctuations close to the surface as that of the 
atomistic model. The use of more detailed potentials is fully com¬ 
patible with the SCE formulation developed. In the inset to Eigure 
Eigure 3 the local volume fraction is depicted. The local volume 
fraction and the structural features obtained from the solution of 
the SCE are the same as the ones by Daoulas et al.'^' Eor more de¬ 
tails and a thorough study of several structural features the reader 
is referred to that work. 

Adhesion Tension 

Within the maximum term approximation, which leads to Eq. ( 6 ) 
for the effective Hamiltonian at the optimum, the grand potential 
for the considered system is estimated as: 

a{V,T,^l)=-kBT\nZ 

= -Po j d^rw'{r)0{r) 

+ Po I dViys(r)0(r) 

-||/dV^.(r) (34) 

Eq Eq. (34), applied to the same volume V of homogeneous poly¬ 
mer in the bulk, which will contain «bulk chains at chemical po¬ 
tential fl and temperature T, yields: flbulk L/t) = —Wbulk/jS = 
— (poV)/(Aj 8 ) since w' (r) = 0 and (r) = 1 in the bulk. Introduc¬ 
ing the subscripts “f” and “fs” to distinguish between free polymer 
surface tension and solid/polymer interfacial tension, respectively. 



_I_I_I_ 

0 0.5 1 1.5 2 


z (nm) 

Figure 3. Self consistent field per segment, W (z) / {Nk^T) for the two 
different chain lengths considered, N = 50 and N = 100. In the inset to the 
figure, the local volume fraction is also depicted. 


and denoting the surface tension of the clean (pure) solid as 7 s, we 
define the adhesion tension^® as the difference 7 s — 7 fs. Now, con¬ 
sidering the special case of a polymer melt next to an undeformable 
solid surface, 

a (Ts - 7fs) = i^buik {V,T,p)-Q. (T, T, p) (35) 

where a is the total interfacial area of contact. The asymptotic 
value obtained from the solution of the SCE problem is ( 7 s — 7fs) = 
73.6 mN/m, which is in good agreement with the adhesion ten¬ 
sion estimated by atomistic MD simulations, ( 7 s — 7fs)lMD = 70± 
10 mN/m 

To the best of our knowledge, experimental data on interfacial 
thermodynamic properties of the polyethylene/graphite system are 
not readily available. According to the theory of Girifalco and 
Good^® the adhesion tension, can be expressed in terms of the ge¬ 
ometric mean of the surface tensions of the liquid and the solid as: 

7s-7fs~2<l>(7s7fs)'/^-7f (36) 

The experimentally measured surface tension of PE at 450 K 
is = 28.1 mN/m,®®’^* and that of graphite is 7 “*’ = 
115 mN/m.®^“®"^ Based on these values, the experimental estimate 
of the adhesion tension, taking the interaction parameter <5 = 1 , is 

(ys-yfs)exp = 85.6mN/m. 

Propagators of End-Constrained Chains 

In order to calculate the configurational integral of end-constrained 
chains based on Eq. (22), we have introduced the propagator 
G>j (z( 0 ),z/^) which is proportional to the probability that a chain 
which has started at height z( 0 ) is found at height z/ after contour 
length J, the whole contour of the chain from 0 to v lying entirely 
in the region of space ^ > z. The propagator is obtained as the 
solution to a diffusion equation ( Eq. (20)) in the presence of the 
W (z) field and is presented in Eigure Figure 4. The chain length 
considered is A = 50, while the end of the chains lies at position 
z(0)=Ng. 

For relatively large values of z', G>j (z(0),z/j) obeys a Gaus¬ 
sian distribution, whose height and breadth increase as z decreases. 
At low z' values, it drops steeply to zero due to the boundary con- 
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dition G>j (z(0),z' = z,s) = 0, yielding an overall asymmetry. For 
the special case z = 0 the propagator (z(0),z',s) resembles a 
Gaussian distribution as long as the chain does not come close to 
the solid surface. Under the influence of the substrate, the probabil¬ 
ity of finding a segment close to it exhibits a local maximum. Thus, 
for z = 0 and intermediate values of s, the distribution appears to be 
bimodal, having one maximum close to z! = z(0) = Rg where the 
z(0) end has been anchored, and a second one inside the attractive 
region close to the substrate. The distribution becomes broader as 
the contour length, 5, increases. This is expected since the chain ex¬ 
plores more space, moving away from its end. For z = 0 the chain 
segments are likely found inside the adsorption layer, contributing 
to the local maximum of the G>j; (z(0),z',^) close to the solid sur¬ 
face. On the contrary, if the dividing surface is placed at higher 
z > 0, the propagator is forced to become zero (due to the adsorb¬ 
ing boundary condition at z' = z), implying also that the probability 
of finding chain segments close to the surface drops to zero. 

Configurational Integral of End-Constrained 
Chains 

Next, we consider the configurational integral ^(z)/^ulk which 
is proportional to the probability that all segments of a chain 
find themselves at a height greater than z from the substrate. 
^ (z) / i^bulk is expected to be a decreasing function of z, since more 
polymer chains are excluded upon increasing z. Moreover it goes 
smoothly to zero at z = min (z(0), z( 1)), since at least one end of the 
chain is found at z. In Figure Figure 5, the dependence of this quan¬ 
tity on the distance of the ends of the strand from the solid substrate 
is presented. We consider two chains of length N = 50 and = 100 
carbon atoms, whose ends are separated by a distance of 0.5Rg par¬ 
allel to the surface. We observe that the profile of ^ (z) exhibits 
a spike at z = 0. This feature indicates strong adsorption of the 
chains onto the solid surface. This effect is more pronounced if the 
self-consistent field w (r) consists only of the wall potential. Us (r). 
If the polymer field is switched off, chains are extremely strongly 
attracted by the wall. This results in extremely an overwhelmingly 
high probability of getting adsorbed and monotonically decreasing 
configurational integrals. At the point where both fields (either SCF 
or wall) approach zero, the results are identical. 

For a given distance from the surface, z, we observe that the prob¬ 
ability of finding a chain for which all segments lie above this dis¬ 
tance increases with the height at which the ends are constrained. 
This is reasonable, since the higher the ends lie from the surface, 
the more difficult it is for the chain to approach the solid surface, 
thus increasing the probability of finding its segments away from it. 
Despite the similarity of the profiles, there exist slight differences 
between the two different chain lengths considered. The differences 
of the profiles vanish as we move the ends of the strand higher. This 
can be attributed to the fact that the range of the attractive potential 
is independent of the length of the chains, thus affecting differently 
the two chain lengths. At this point, we have to stress that the con¬ 
figurational integrals (as presented in Rg units in Figure Figure 5) 
are completely different for the two chain lengths, if considered in 
the real spatial scale of the problem. 

Equilibrium Probabilities 

We now consider the equilibrium probability, P^, of finding a chain 
in the adsorbed state, as this was defined in terms of the configura¬ 
tional integral ( Eq. (24)). In Figure Figure 6 the contour plots of 
the probability to find a chain in the adsorbed state are presented 
for the two chain lengths considered, N = 50, (a) and N = 100, 



0 0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 2 

z (flg units) 

Figure 5. Configurational integral of end-constrained chains, 3f{z)/ f^bulk; 
as a function of distance from the solid surface, for two chain lengths. The 
ends of the chains are progressively elevated from the solid substrate, start¬ 
ing from z(0) = z(l) = 0.5Sg to z(0) = z(l) = 2.0Sg. The distance between 
the ends of the strand at an ^cy-level parallel to the surface is O.SRg. 


(b), respectively. The similarity of the configurational integrals is 
further manifested in the adsorbed state probability. The contours 
are almost identical, given that the spatial scale is expressed in Rg 
units. For comparison, the real distances (in nm) are also noted 
in the plots. For distances of the ends from the surface smaller 
than or close to Rg, the probability of finding the chain in the ad¬ 
sorbed state (i.e. with at least one of its segments lying below 
z = 5) is almost unity, due to the proximity to the attractive sur¬ 
face. For larger distances, the probability of finding the chain in the 
adsorbed state is reduced, exhibiting a roughly hyperbolic depen¬ 
dence on the distances of the ends of the strand from the surface. 
For z(0),z(l) — 2Rg it is highly unlikely to find an adsorbed chain. 
If one of the two ends of the chain remains at a short distance from 
the surface (e.g. up to Rg), there is always a finite probability of 
finding the chain in the adsorbed state, irrespectively of the posi¬ 
tion of its other end. This is a consequence of modeling chains 
as infinitely extensible Gaussian threads. Finally, the contour plots 
are symmetrical about z(0) = z( 1). This stems from the fact that the 
choice of which end of the chain corresponds to v = 0 and which 
corresponds to v = 1 is arbitrary. 

Adsorption Potential of Mean Eorce 

Based on the configurational integral, J3f(z) /^ulk> we have defined 
a free energy per unit length, A(z), of a chain to have at least one of 
its segments below a certain distance, z, from the substrate. To com¬ 
pute A(z) we apply a spline interpolation to JSf(z) in order to obtain 
a smooth estimate of its first derivative, d.^(z)/dz, which we then 
use in order to estimate the free energy profile, A (z), via . Eq. (26). 
The resulting profiles for N = 50 and various distances of the ends 
of the chain from the substrate are depicted in Eigure Figure 7. The 
normalization constant C 4 which is needed in Eq. (26) has been set 
to 1 m^ in order to cohere with the units of m^ in which the configu¬ 
rational integral has been expressed. The potential of mean force is 
expected to become infinite for z min (z( 0 ), z( 1 )), due to the hard 
constraint we impose with respect to the height of the ends from the 
surface. This behavior is evident in Eigure Figure 7, where we have 
anchored both ends at the same height, for simplicity. 

The free energy profiles exhibit a well defined minimum ex- 
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z' (R„ units) 



z' (Rg units) 
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Figure 4. The propagator G>z (z(0),z',i) of a chain of length = 50 to find its whole contour from z(0) to i above the z level, as a function of z'■ Four 
different values of the height z have been considered, which are plotted in different subfigures. For a given z, the distribution becomes broader as i increases. 
The end of the chains has been anchored at height z(0) = Rg. Please note the different scale of the ordinate of the subfigures. 



Figure 7. Potential of mean force, A(z), for an end-constrained strand 
to have its lowest bead at distance z to z -f dz from the solid substrate ( 
Eq. (26)). The chain consists of 50 CH 2 groups and its unperturbed radius 
of gyration is Rg = 1.27 nm. The normalization constant used in Eq. (26) 
was C 4 = 1 m'*. In the inset to the figure, the potential of mean force in the 
region 0 < z < 0.4 nm is displayed in order to facilitate the estimation of the 
depth of the minimum of A(z) close to the wall. 


tremely close to the solid surface (z = 0), whose depth depends 
on the distance of the ends of the chain from the surface. For 
z(0) = z(l) < l.ORg the adsorption of the chains on the surface is 
spontaneous, the free energy with respect to distance being a mono- 
tonically increasing function. The only thermodynamically stable 
state for the chains is the adsorbed one, in agreement with the ad¬ 
sorption probabilities presented in Figure Figure 6. The fact that no 
free state exists for these distances from the wall is reminiscent of 
the behavior of terminally grafted chains in ref 42, which are shown 
to be exclusively in an adsorbed state when the segment-substrate 
attractive interaction energy exceeds a certain critical value. For 
distances of the ends z(0),z(l) > Rg, a second minimum of the 
free energy appears, corresponding to a “free” or desorbed state of 
this chain in the bulk. Thus, a barrier separating the two minima 
appears. The height of the barrier grows as the ends of the chain 
are anchored further apart from the solid surface. This is expected, 
since it becomes progressively more difficult for the chains to ex¬ 
pose their segments close to the surface. Moreover, the position 
of the minimum shifts to larger distances, indicating that the ad¬ 
sorption becomes even more difficult. On the contrary, the position 
of the barrier seems to remain unaffected upon increasing the dis¬ 
tance of the ends from the solid surface. This facilitates the clear 
definition of a dividing surface centered at the top of the barrier, 
8 ~ 0.3 nm, relative to which adsorption and desorption are con¬ 
sidered to take place. The maximum of the free energy coincides 
with the width of the square-well potential (Figure Figure 3) and 
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Figure 6. Probability of finding a chain in the adsorbed state (at least one segment in the region z < 5) as a function of the distance of the ends of the chain 
from the graphite surface. The end-constrained strand consists of either 50 (a) or 100 (b) methylenes. The position of the dividing surface, 5, has been obtained 
as the maximum of the A(z) profiles (see Figure Figure 7). 


is independent of the length of the chains. The profile at the top of 
the barrier can be locally approximated by a polynomial function of 
high degree (6 or 8 ), allowing us to precisely define its maximum 
(dA(z)/dz|,^g = 0 ) and obtain an accurate estimate of its curva¬ 
ture (d^A(z)/dz^|_g) needed in eqs Eq. (27) and Eq. (28) in order 
to estimate the adsorption and desorption rate constants from eqs 
Eq. (30) and Eq. (31), respectively. 

Adsorption and Desorption Rates from SCF/TST 

The last and most crucial step in our study is the estimation of 
adsorption and desorption rate constants for an end-constrained 
strand. The identification of the free energy barrier height (Eig- 
ure Eigure 7) and the estimation of its curvature allows us to ap¬ 
ply dynamically-corrected transition state theory to extract the rate 
constants for a strand to overcome the barrier of free energy, in 
order to either adsorb or desorb. In the following, we will use 
Cstrand = ('”strand/™CH2 ) ?0 the Mction coefficient of the strand, 
invoked in eqs Eq. (30) and Eq. (31). At first, in Figure Figure 8 
we restrict ourselves to equidistant ends of the strand from the solid 
surface. The rate constants we present are independent of the dis¬ 
tance between the projections of the ends on the xy plane, due to 
the factorizability of the propagators ( Eq. (19)) that leads to can¬ 
cellation between numerator and denominator of the last terms in 
eqs Eq. (30) and Eq. (31). Based on Eq. (30) we estimate the ad¬ 
sorption rate constant of chains of length 50 and 100. The results 
are presented together with the estimates of adsorption rate constant 
obtained by hazard-plot analysis of MD trajectories for the same 
system (see next subsection). 

The adsorption rate, as presented in Figure Figure 8 , exhibits an 
extremely strong dependence on the distance of the strand’s ends 
from the solid surface. The nearly exponential decay of the rate 
constant with distance of the constrained ends from the surface can 
be observed for both chain lengths under consideration. In the inset 
to Figure Figure 8 , the same data are presented with the abscissa 
measured in Rg units. For both chain lengths, the estimation of 
the rates can be accomplished down to 0.57?g; below this lower 
limit adsorption is spontaneous, leading to infinitely high adsorp¬ 
tion rates. At intermediate distances, a nearly exponential decay of 



Figure 8. Adsorption rate constant, as a function of the distance of 
the ends from the graphite surface. Both ends are kept at the same dis¬ 
tance from the graphite surface, z(0) = z(l) and we consider chain lengths 
of A = 50 and N = 100 methylenes. In the inset to the figure the same 
data are presented with the abscissa in Rg units. Adsorption rates for an 
end-constrained Cjq subchain in a C 200 polyethylene melt next to graphite 
obtained by MD simulation using a detailed united-atom model are also 
shown (compare Figure Figure 11). 
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the adsorption rate constant takes place, for both chain lengths con¬ 
sidered. At even larger distances, a downward departure from the 
exponential trend can be observed. The adsorption rate should go to 
zero when the triangle formed by the end-to-end distance and two 
fully extended substrands of total length equal to the chain length, 
forming equal angles with the surface normal, barely reaches the 
surface. The infinite extensibility of the Gaussian model may ac¬ 
count for the gentler drop of the adsorption rate constant, as ob¬ 
served in Figure Figure 8. The adsorption rate constants depend 
also on the molar mass of the strands. At small distances from the 
surface, shorter, thus more mobile, chains can approach the sur¬ 
face more easily. However, at larger distances, where stretching of 
the chain should take place to reach the surface, longer chains with 
more segments adsorb faster. The time scale at which adsorption 
takes place cannot be fully accessed by MD simulations, especially 
for the longer chains. However, where results are available from 
direct measurement of adsorption rates, the agreement between the 
SCF/TST approach and the MD simulations is favorable. 

Rate constants for adsorption/desorption from the melt are dif¬ 
ficult to measure directly. Dietsche et al.^^ have measured the 
displacement of deuterated alkanes by hydrogenous alkanes in an 
C[g alkane liquid flowing past a surface (ZnSe, Cr, Cu) with at¬ 
tenuated total reflectance Fourier transform infrared spectroscopy 
(ATR/FTIR). The ATR/FTIR experiments revealed a surprising 
interaction between the weakly adsorbing oligomeric polyolefins 
and the surface, which manifests itself in slow exchange between 
the surface and the bulk. The picture which seems to emerge is 
an entropically-driven flow-rate- (or stress-) dependent adsorption 
equilibrium which retards the exchange between the surface and 
bulk. In order to reproduce the experimental findings with finite- 
element calculations, these authors have assumed a diffusivity of 
3.1 X 10^*® m^ s^*. If we imagine a strand diffusing in the inter¬ 
facial region with this diffusivity, the timescale of the diffusion of 
a strand over a lenght commensurate with its radius of gyration is 
~ 10^* s, which is relevant to the results obtained by the SCF/TST 
approach. 

The next step. Figure Figure 9, is to allow the distances of both 
ends of the chain from the substrate to vary independently. The 
adsorption rate contour plot. Figure Figure 9, is symmetrical about 
z(0) = z(l), since the choice of which end of the chain corresponds 
to i = 0 and which corresponds to i = 1 is arbitrary. It is evident 
that the rate constants extend over a wide range of timescales (from 
~ lO'^* s^* at 1 nm to ~ 10^ s^* at 5 nm). The adsorp¬ 
tion rate for both chain lengths (Figure Figure 9 (a) and (b)) drops 
quickly when moving away from the surface. It should be noted 
that the adsorption rate drops almost by one order of magnitude 
upon elevating the ends of the strand by 1 nm from the graphite 
surface. However, if one of the ends of the chain finds itself close 
to the solid substrate, the adsorption rate constant remains relatively 
high, irrespectively of the height of the other end. Overall, at the 
same values of z(0) /Rg andz(l)/Rg, longer chains (Figure Figure 9 
(b)) adsorb more difficultly on the solid surface, than shorter chains 
do (Figure Figure 9 (a)). This may be attributed to the increased 
friction they experience due to interactions with their environment, 
resulting in reduced mobility. 

The dependence of the desorption rate on the distance of the 
ends from the substrate is presented in Figure Figure 10. As ex¬ 
pected, desorption rates are increasing functions of the distance 
of the ends of the chain from the surface and exhibit a symmetry 
about z(0) = z(l), since the assignment of the chain ends to i = 0 
and ^ = 1 is arbitrary. Close to the surface the desorption rates are 
lower than the corresponding adsorption rates. As one moves fur¬ 
ther away from the substrate, however, the balance of the two rate 
constants radically changes. Shorter chains exhibit the tendency to 


desorb faster, under the same conditions. For both adsorption and 
desorption rates, the detailed balance condition holds. This means 
that the probability of being free multiplied by the adsorption rate 
constant equals the probability of being adsorbed multiplied by the 
desorption rate constant, for all cases considered. 


Adsorption Rates from MD 


In order to evaluate the rate constants for the adsorption/desorption 
events, we invoke a hazard plot analysisof the atomistic MD 
trajectory. Hazard plot analysis is based on the evaluation of the 
cumulative hazard. The hazard rate, h{t), is defined such that h{t)At 
is the probability that a system which has survived a time f in a 
certain state since its last transition, will undergo a transition at 
the time between t and dt. The cumulative hazard is defined as 
H{t) = jQh{t')dt'. By assuming a Poisson process consisting of 
elementary transitions with first order kinetics from one state to the 
other, the Poisson rate can be extracted as the slope of a plot of the 
cumulative hazard versus the residence time. Thus, the adsorption 
rate constant can be obtained from the hazard plot of the residence 
time of the system in the free state, and vice versa. The dividing 
surface of the adsorbed and the free state, is placed at 5 = 0.3 nm, 
as estimated from the SCF calculations. 



Residence time in the free state, t (ns) 


Figure 11. Hazard plot of the time between transitions from the free to the 
adsorbed state, as calculated atomistically with molecular dynamics. The 
height at which the ends of the chain are anchored is varied. Linear fits are 
also included which provide an estimate of the slope, i.e. the adsorption rate 
constant. 


In Figure Figure 11 we present the adsorption hazard plot for an 
end-constrained intra-chain strand of 50 CH 2 united atoms (such as 
the red one in Figure Figure 2), whose ends have been anchored 
at z(0) = z(l) = 0.7,0.8 and 1.3 nm. H{t) becomes linear at long 
times. At short times it exhibits some nonlinearity, especially when 
z(0) = z(l) is large, betraying fast recrossing events at the dividing 
surface. The slope of the curves corresponds to the rate constants 
for this strand to make a transition from the free to the adsorbed 
state. It is evident that the time a strand spends in the free state 
before adsorbing onto the solid substrate grows quasi exponentially 
with the distance of the ends from the surface. With our compu¬ 
tational resources, a study for anchoring distances larger than 1.3 
nm was not feasible, due to the maximum time that could be ac¬ 
cessed with an MD simulation. Even for smaller distances, there is 
a significant error related to fitting the long time part of curves to 
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Figure 9. Adsorption rate constant, as a function of the distance of the ends of the chain from the graphite suface. The chain consists of either 50 (a), 
or 100 (b) methylenes. The radii of gyration ai‘e /?g(50) = 1.27 nm and i?g(100) = 1.80 nm, respectively. 



A/= 100 (1400 g/mol) 



0.6 1 1.4 1.8 2.2 2.6 3 


(1.44 nm) z(0) (Rg UnitS) 

(b) 


Figure 10. Desorption rate constant, as a function of the distance of the ends of the chain from the graphite suface. The chain consists of either 50 (a), 
or 100 methylenes. The radii of gyration are /?g(50) = 1.27 nm and i?g{100) = 1.80 nm, respectively. 
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a straight line. This is due to the fact that at short time-scales a lot 
of short-lived events are monitored because of the recrossing of the 
dividing surface. The curves depicted in Figure Figure 11 provided 
us with the estimates of adsorption rate constants included in Fig¬ 
ure Figure 8, where the comparison against the SCF/TST estimates 
was made. 

Summary and Conclusions 

In this paper we have outlined a strategy that combines field the¬ 
oretic modeling with transition state theory to compute rate con¬ 
stants for adsorption and desorption at a melt/solid interface of 
a (sub)chain whose ends are constrained in space, e.g. through 
crosslinks or entanglements. At first, the problem of determining 
the equilibrium properties of a polymer melt of specific chemical 
constitution, adsorbed on a certain solid substrate through a contin¬ 
uum SCF approach was considered, following the work of Daoulas 
et al. The Gaussian chain model reproducing polyethylene chains 
has been used. As far as the interactions with the substrate are 
concerned, these were taken into account via a square well po¬ 
tential whose parameters were matched to reproduce the graphite- 
polyethylene interaction. The interfacial structure and thermody¬ 
namics were obtained from the solution of the SCF problem. An 
estimate of the adhesion tension is provided, for the first time to the 
best of our knowledge, which has been found to be in reasonable 
agreement with that computed by a realistic united-atom simulation 
and with available experimental evidence. 

Based on the solution of the SCF problem, a conveniently de¬ 
fined propagation function allowed us to estimate the configura¬ 
tional integral of end-constrained chains with a given distance of 
closest approach to the solid substrate. From the configurational 
integral we computed the Helmholtz energy (potential of mean 
force) of the end-constrained chain as a function of our chosen one¬ 
dimensional order parameter (the distance of closest approach to 
the solid surface). Chains having their ends constrained close to the 
substrate exhibit a monotonically increasing free energy when mov¬ 
ing away from it, with a well-defined minimum close to it. Thus, 
spontaneous adsorption of these chains takes place. Chains whose 
ends are constrained at larger distances from the substrate exhibit 
a free energy profile with two well-defined minima (states). The 
free energy barrier in between the two minima serves for defining 
a dividing surface which discriminates between “adsorbed” (close 
to the wall) and “free” (far from the wall) states. The distance of 
the barrier from the substrate surface is very close to the width of 
the square well potential used for representing the graphite/polymer 
segment interaction and is independent of the chain length. 

Dynamically corrected transition state theory, in the form of 
Kramers theory in the high friction limit, applied to the free energy 
profile, provides a reasonable way of describing the dynamics of 
adsorption/desorption of the end-constrained chains. The friction 
factor needed for calculating rates is connected to the monomeric 
friction coeffecient, which has been extracted from atomistic MD 
simulations in the Rouse regime.®^ The adsorption/desorption rate 
constants are found to depend strongly on the distance of the ends 
of the strands from the solid surface. The length scale of adsorption 
is limited to a few nanometers for chains composed of 50 to 100 
methylene units. Despite the locality of the effect, a decay by four 
orders or magnitude (or more) of the adsorption rate constant is ob¬ 
served when moving from to 1 to 5 nm from the graphite surface. 
The rate constant for adsorption decays practically exponentially 
with the distance of the chain ends from the solid surface over an 
intermediate range of distances between the longest distance where 
spontaneous adsorption takes place and three times the unperturbed 
mean square radius of gyration. For larger distances of the ends 


from the surface the rate constant for adsorption drops more pre¬ 
cipitously. The rate constant for desorption increases with distance 
between the chain ends and the solid substrate. It is generally less 
sensitive to distance than the adsorption rate constant. 

To the best of our knowledge, this the first time that the solution 
of the polymer SCF problem next to a solid surface has been used 
to provide insight into the dynamics of the system. This work sug¬ 
gests that SCF, coupled with transition-state theory, can give good 
results for the rate of adsorption/desorption between polymers in a 
melt and a surface. SCF theory is less computationally demanding 
than atomistic simulations, and thus SCF can be used to examine 
regimes not practically accessible to atomistic simulations, such as 
the regime of low adsorption rates. Furthermore, the parameters in 
SCF have been directly connected to atomistic parameters, suggest¬ 
ing that SCF may provide a more faithful description of the problem 
than alternative, coarse-grained simulation methods. In the regime 
where molecular dynamics simulation of adsorption rates is prac¬ 
tical (short chains, very close to the surface), our results for the 
adsorption rate are in good agreement with the results from molec¬ 
ular dynamics simulation. Finally, a reasonable agreement between 
the melt adhesion tension computed from SCF and that measured 
in atomistic simulations and in experiment, was found. 

In the ongoing debate of interfacial dynamics, our results suggest 
that the crosslinks or entanglements that are nearest to the surface 
would be the controlling factors influencing the surface dynamics. 
Some of these links may be permanently attached to the surface, 
while other may adsorb/desorb with a rate constant. If the magni¬ 
tude of the adsorption rate constant is larger than the time needed to 
entagle/disentagle, the surface may enhance the entanglement den¬ 
sity of the polymer, thus contributing to its mechanical properties. 
If the two surfaces come closer, in a way that segments lying in the 
middle can adsorb on both of them, bridges of adsorbed may be 
formed, slowing down the dynamics of the polymer in the slit. 
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Appendix A: Grand Partition Function of an Ensemble of Noninteracting Chains Subject to a Field 


We consider a random fluctuating field w (r) and form the functional integral: 


^[j3w(r)]exp| I dV -i (0 (r) - l) jSpow( 




with i being the imaginary unit, j3 = 1/ (ksT) and the volume fraction operator 0 (r) given by Eq. (2). We can rearrange the internal integral 
and write the functional integral as: 


/^[/3w(r)]exp|-/d3r[(^) ' Po-W + 


exp — / ^-^d r 


The functional integral over w (r) in the last expression is Gaussian and should yield a number C ^ as a result, which is dependent upon po, 
k:j, and T. Thus, we can write 


exp ^J ^^2KlkBT =C I Si[Pw{r)]expi^Jd?r -i (r) - l)/3po 


(^Pow(r))3 


Substituting the left-hand side of Eq. (A.3) with its right-hand side within Eq. (1) we obtain: 


= Ey|exp(0)iV«c|l?[i3w(r)] f n l?r„ (.) ^ [r„ (.)] 
expj /d^r -i/3poH'(r) (^(r)- l) - ^l/s (r) ^ (r) - (/3pow(r))3 


and, recalling Eq. (2), 


S = £^-p(^0)A«c/l?[/3>v(r)] 

n |/ ^rcf (•) .^ [r« (•)] exp - J <^r As 8 {r-ra{s)){PiNw{r) + fiNUAr)) 


expj/d^r i/3pow(r)-^^5:^(/3pow(r))3 | (A.5) 

In Eq. (A.5), the integration over (•) is reduced to a partition function of n independent Gaussian chains, the segments of each chain 
interacting with the field !w(r) + 17s (r)- Introduction of the fluctuating field H’(r) has served to decouple the chains and reduce the problem 

of n interacting chains to a problem of n single chains interacting with an effective field in the calculation of the grand partition function. 

Eollowing Edwards, we introduce the notation Q[iw + Us\ to indicate the partition function of a single chain subject to the field iw + Us, 
acting on its segments, relative to the partition function of a field-free chain. Then, 

J ^r„(»),^[rc((»)]exp d^r^ d^ 5 (r-(j)) (jSi Aw(r) + j3A(/s (r)) 

= |y {•)^[ra (•)]|2[(>V-l-(/sl 

=Zf,,,Q[iw + U,] (A.6) 

where Zf^e = f Sir a (•) [ra (•)] is the partition function of a free chain. With this notation, 

E = c£ yexpT^^y (iVZfree)" /^[/3H’](e[itV-f f/s])" 

expjyd^r ij3pow(r)-(j3pow(r))3 | (A.7) 


: j S[fiw\expi^- j A^r l5pow{r) +{ppow{r)f j; 

£ (^2free)"(e[!W-b(7s])" 


Now, the summation over different n can be performed, yielding 


£ ^exp^^^y {NZf,ee)"{Q[iw + Us])'‘ 
n=o”- / 


= exp|exp y — jNZf,ssQ[iw+Us] 


which leads to the expression of Eq. (4) for the grand partition function. 
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Appendix B: Saddle-point Approximation 


We replace the functional integral appearing in Eq. (4) of the main text for E with its dominant term, obtained by setting the functional 
derivative of H to zero, SH/8w = 0: 


^ = ^ipo + tCTPow(r) - ^ exp 


Now, by definition,®^ Eq. (A.6) of Appendix A can be rewritten as: 


Q[iw + Us] = 


I ^rc((»),^[r„(»)]exp|-y^d^r^ d.s 5 (r-r„ (j)) (/3f Aw (r)+/3Af7s (r)) 
f&ra{»)i^[ra{»)] 


From Eq. (B.2), 


SQ^ _ 

/^r„( 


- / 

Jd?rJ d.v5(r 


(•) .^ [r„ (•)] -jSiA / di-5(r-r„(5)) x 


Lv 5 (r - r„ (j)) (j8f Aw (r) + fiNUs (r)) 


= -jSiAg< &ra{») 3^[ra{»)] / d.?5(r-ra(j)) x 


exp — / d^r / d^S (r — Ta (i)) (j8iAw(r) + PNUs (r)) 


/|y^rc((*).^[ra(»)]exp l^-y dVy di5(r-ra ( 5 )) (/3iAw(r) +/3Af7s (r)) | (B.3) 

Recalling the definition of ^ from Eq. (2), Eq. (B.3) can be written as: 

ib.4) 

with () denoting an average over the distribution defined by fhe single-chain partition function Q[iw + Us], Eq. (B.2). In averaging 0 (r) of 
Eq. (2), all n chains involved in the summation for a = 1 to m yield the same contribution. 

We introduce the symbol 

, z s / ^ z ^ _ A _ 


r /*! 
J Jo 


0(r) = (0(r))^^-= — dsd{r-ra{s))j 

with n being the total number of chains present in the considered interfacial region at the optimum. Then, Eq. (B.3) gives: 


ow n 


Substituting Eq. (B.6) in Eq. (B.l), we obtain: 


^ =-'Po + K'TPo'v(r) + f^exp("^)AZfteee[iw + f/s](|)(r) =0 


n \kBT 


whence the self-consistent field is obtained as: 


w (r) = —^ 11 - t exp r ^']NZfyeeQ [iw + (/s] 0 (r) 
'^tPo I n yksTJ 


Very far from any interfaces, bulk conditions will prevail: 


w(r) = 0 
0 (r) = l 


By applying Eq. (B.8) in that region, we are led to the condition: 


fi \kBTj 

Using Eq. (B.IO) in Eq. (B.8), we obtain the self-consistent field as 


^exp { ^)AZfreee[fw-|-f7s] = 1 


w (r) = -—[!-()) (r)] 
KtPo 


At the optimum. 


dVpow(r) = ^ dV[l-(|)(r)] = ^ ( V - ^ ) = ^ 


i nN\ i 


J KtJ Kt \ Pq j Kt PO 

with Mbulk the number of chains within volume V of the bulk polymer under the considered conditions of T and pA. 
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FromEq. (B.ll), 

[pow(r)]2 = - ^ Jd^r[^ (r) - 1]^ (B.13) 

In view of eqs Eq. (B.12) and Eq. (B.13), the first sum in Eq. (5) of the main text for the field-dependent effective Hamiltonian at the optimum 
becomes: 

ydV[-ipoiv(r) + y (poiv(r))^] = ("bulk-") - (r) - 1]^ (B.14) 

Following a different route, using eqs Eq. (B.ll) and Eq. (B.13), and the fact that po = "bulk^/^ we obtain an alternative expression for the 
same quantity as follows: 

ydV[-ipow(r) + y (pow(r))^] =_"^^ jd?riw{r)^{T) 

+ ^ydV[l-(/.(r)]2 (B.15) 

By virtue of Eq. (B.15), the effective Hamiltonian, Eq. (5) can be written as: 


where, from Eq. (B.IO), 


And by setting 


jd riw{r)(j>{r) + —J d r[l-(j>{r) 


^exp (^^^NZf,ssQ[iw + [/s] = j 


w' (r) = i w (r) -f Us (r) 


which is a real field, we arrive at the effective Hamiltonian as expressed in Eq. (6) of the main text. 


(B.16) 

(B.17) 

(B18) 
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